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q " Abstract 

We investigate the history of quantum measurements on scattered 
| probes. Before scattering, the probes are independent, but they be- 

come entangled afterwards, due to the interaction with the scatterer. 
The collection of measurement results (the history) is a stochastic pro- 
cess of dependent random variables. We link the asymptotic properties 
of this process to spectral characteristics of the dynamics. We show 
that the process has decaying time correlations and that a zero-one law 
holds. We deduce that if the incoming probes are not sharply localized 
with respect to the spectrum of the measurement operator, then the 
process does not converge. Nevertheless, the scattering modifies the 
measurement outcome frequencies. We show that those are the fluxes 
l/-) | of the von Neumann measurement projections in the initial probe- 

scatterer state. We illustrate the results on the Jaynes-Cummings 
^ . model. 

d 

1 Introduction and main results 

We consider a scattering experiment in which a beam of probes is directed at 
a scatterer. The probes are sent to interact sequentially, one by one. Before 
the scattering process, they are identical and independent. The interaction 
of each probe with the system is governed by a fixed interaction time r > 
and a fixed interaction operator V. After interacting with the scatterer, 
a quantum measurement is performed on each "outcoming" probe. The 
result of the measurement of the n-th probe is a random variable, denoted 



03 



> 



*Department of Mathematics and Statistics, Memorial Univesity of Newfoundland, St. 
John's, NL Canada, A1C 5S7; merkli@mun.ca; Supported by an NSERC Discovery Grant 

t Current address: Department of Mathematics, University of Toronto, Bahen Cen- 
tre, 40 St. George St., Toronto, ON Canada, M5S 2E4; mark.penney@mail.utoronto.ca; 
Partially supported by NSERC USRAs 



X n . The stochastic process {X n } n >i is the measurement history. Due to 
entanglement of the probes with the scatterer, the X n are not independent 
random variables. We analyze asymptotic properties of this process. 

A concrete physical setup is given by atoms (being the probes) shot 
through a cavity containing an electromagnetic field, the modes which in- 
teract with the atoms forming the scatterer. We assume that the incoming 
probe states are stationary with respect to their isolated dynamics. 

We study systems with only finitely many degrees of freedom involved 
in the scattering process. This means that the Hilbert spaces of pure states 
both of the system and each probe is finite-dimensional. The measurement of 
a probe is a von Neumann, or projective, measurement associated to a self- 
adjoint probe measurement observable M. The eigenvalues m of M are the 
possible measurement outcomes. Due to finite dimensionality, the random 
variables X n have finite range. 

The present work can be viewed as the continuation of recently developed 
techniques for the mathematical analysis of repeated interaction quantum 
systems [12, 13, 14, 15]. In these references, asymptotic properties of the 
scatterer have been investigated, without considering the fate of the outcom- 
ing probes, and without quantum measurements. While the setup of our 
present work is similar to the one in the given references, our focus here is on 
the measurement outcomes process. We show that generically, this process 
does not converge. We describe the fluctuations on the measurement history, 
provoked by the scattering process, by analyzing the measurement frequen- 
cies. A more detailed comparison to related works is given at the end of this 
section. 

As explained in the references above, in absence of quantum measure- 
ments on probes, and under a generic ergodicity assumption, one shows that 
the scatterer approaches a so-called repeated interaction asymptotic state af- 
ter many interactions. We keep this assumption in the present work. 

(A) Assume that if no measurement is performed (M = 11), then, under the 
repeated interaction with the probes, the scatterer approaches a final 
state. The convergence is exponentially quick in time. 

The precise mathematical formulation of this assumption is given in Sec- 
tion 2.4, see before (25). It is a condition on the spectrum of a reduced 
dynamics operator, and necessitates the introduction of some technicalities 
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which we want to avoid in this introduction. Condition (A) is generically 
satisfied, and is not hard to be verified explicitly, and one even calculates the 
rate of convergence for concrete models (see the above references). 

We now explain our main results. Denote by cr(X r , . . . ,X S ) the sigma- 
algebra generated by the random variables X r , . . . ,X S , 1 < r < s < oo. We 
denote by P the probability measure associated with the process {X n } n >i. 

Theorem 1 (Decay of correlations) Suppose that Condition (A) holds. 
There are constants c, 7' > 0, such that for \ <k<l<m<n< 00, 
A G cr(Xk, ■ ■ ■ Xi) and B G o-(X m , . . . , X n ), we have 

\P{A nB)- P(A)P(B)\ < cP{A) e -' 1 ' {m ' l) . (1) 

We give a proof of the theorem in Section 4.1. Intuitively, the system starts 
relaxation to its asymptotic state during the time m — I between two consec- 
utive measurements, and hence erases correlations bewteen the two measure- 
ments. The rate 7' in (1) is linked to the convergence rate in Assumption 
(A), see Section 2.4. 

The tail sigma- algebra is defined by T = n n >ia(X n , X n+1 , . . .). Decaying 
correlations imply the following zero-one law. 

Corollary 2 (Zero-one law) Assume that Condition (A) holds. Any tail 
event AeT satisfies P(A) = or P(A) = I. 

In textbooks, the Kolmogorov zero-one law is usually presented for inde- 
pendent random variables [8]. However, an adaptation of the proof yields 
the result for random variables with decaying correlations, see [1] (and also 
Section 2.4). The tail sigma- algebra captures convergence properties. For 
instance, given any outcome m G spec(M), the set {lim n X„ = m} is a tail 
event, hence, according to Corollary 2, it has probability zero or one. 

We now explain why fluctuations in the process persist generically, for all 
times. Let u in be the state of the incoming probes, denote by E s the spectral 
projection of the measurement operator M associated to S C spec(M) and 
denote E m = E{ m y for m G spec(M). In absence of interaction (V — or 
t = 0), the Xj are independent random variables. We show in Proposition 12 
that the dependence generated by the interaction with the scatterer is small 
for small interactions, uniformly in time. Therefore, since P(X n = m) = 



3 



ou in (E m ) +0(\\V\\), we have P(X n+1 =m,X n = m) = P(X n+1 = m)P{X n = 
m) + 0(|| V ||), and consequently, 



The numbers u)- m (E m ) are probabilities. Thus, ^ m uf n (E m ) = 1 if and only 
if for a single m we have oj- ni (E mo ) = 1 while for all other m, u- in (E m ) = 0. 
This means that P(X n+1 = X n ) < 1 for small V, whenever there are sev- 
eral m with oj\ n (E m ) > 0. Together with the zero-one law, this implies that 
P(X n converges) = whenever the incoming state is not localized in a sin- 
gle subspace of M (and V is small enough). If m is a simple eigenvalue 
of M with associated eigenvector ip m , then w in (E m ) = 1 is equivalent to 
^in(-) = {ipm, -ipm)- Statistical fluctuations in the incoming probes (mixture 
of states localized w.r.t. measurement values) thus get transferred to out- 
coming probes, even in the limit of large times. The following is a more 
general statement of this fact. 

Theorem 3 Assume Condition (A) holds. There is a constant C s.t., for 
any S C spec(M) with u- m (E s ) ^ 1, if \\V\\ < C(l - u in (E s )), then P(X n e 
S eventually) = 0. 

The result on non- convergence of X n explained before Theorem 3 is a special 
case of Theorem 3, when S = {m}, m G spec(M). We mention that our 
analysis also gives a condition under which P(X n e S eventually) = 1, see 
Lemma 11. 

The process X n carries information about the scattering process, encoded 
in the relative occurrence of a particular measurement outcome. We define 
the frequency of m G spec(M) by 



f m is a random variable and the limit (whose existence will be shown) is in 
the almost everywhere sense. The following result analyzes the influence of 
the scattering process on the frequencies. It shows that the main contribution 
is expressed as the flux of the observable 
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f m = lim —{number of k e {1, . . . , n) s.t. X^ = m}. 



n— >oo fl 
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in the initial state us ®uJm- E m {r) is the average of the operator E m over an 
interaction period r. 

Theorem 4 (Frequencies) Assume condition (A) holds, and suppose the 
interaction operator is of the form XV , where A G K. is a coupling constant. 
Then f m = f m (X) exists as an almost everywhere limit and is deterministic 
(not random). It is a holomorphic function in A at the origin, and 

f m {X)=u 1Il {E m ) + Xf' m + 0{X 2 ), 

where f' m = u s ® w in (ir[\/,E m (r)]) . 

Note that by invariance of ug ® oo- in under the free dynamics generated by 
H = H s + H P , we have u s <S> u in (iAr[V,£' m (r)]) = tu s <S> w in (i[H,E m (r)]) , 
where H = H + XV. Thus 

at t=o 

which is (r times) the change in the quantity E m {r) per unit time, i.e., the 
flux. 

Theorem 5 (Mean) LetX n = -(Xi + • • • + X n ) be the empirical average 
of the process {X n }. We have a law of large numbers, 

limX n = /ioo, 

n— s-oo 

almost everywhere, for some deterministic (non random) fi^ G M. // the 
interaction is of the form XV , where A G R is a coupling constant, then 
f'oo = f'oo(X) is a holomorphic function in A at the origin, and 

/ioo(A) = a; in (M) + A /U / 00 + 0(A 2 ), 

where fj,'^ = us <S> Wi n (ir[V, M (r)]) «s £/ie o/ t/ie averaged measurement 
operator M(r) = £ /J" e is ^ Me~ isHv ds . 



Relation to previous work. The literature on repeated interaction 
systems we are aware of can be classified into two categories. In a first one, ef- 
fective evolution equations are derived by taking continuous interaction limits 
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[6, 7, 3, 4, 30, 31, 21] and in a second category, the dynamics is left discrete, 
but the time-asymptotics is investigated [12, 13, 14, 15, 9, 10, 22, 28]. To our 
knowledge, repeated interaction systems have been first proposed in [6, 7], 
as approximations for system-environment type models. It is proven that 
the discrete evolution converges to that of a quantum Langevin equation, in 
the limit of ever shorter system-probe interaction times. The study of the 
time-asymptotics was initiated in [12] (and further developed in [13, 14, 15]), 
where the goal was not to reach a continuous-time limit, but rather to deter- 
mine the time-asymptotics of the system state. It is shown that the (reduced) 
system state converges to a final state. The latter is described explicitly (by 
rigorous perturbation theory in the system-environment coupling) and its 
physical (thermodynamic) properties are found. The present work falls into 
the second category. It would be interesting to examine how our mathemat- 
ical approach leads to effective evolution equations, but we have not studied 
this yet. 

The mathematical formalism we have developed in [12, 13, 14, 15] is 
based on a spectral approach to the time-asymptotics of open quantum 
systems (algebraic formalism of quantum mechanics, Liouville generator of 
dynamics, spectral analysis of the latter - via spectral deformation in the 
infinite-dimensional situation). The present work extends this formalism to 
the setting of repeated probe quantum measurements. The references [9, 10] 
also deal with such measurements. Their analysis is based on the theory of 
classical stochastic processes (law of large numbers, martingale convergence, 
large deviation principle). The two mathematical approaches are completely 
different. The models are different too, and so are the questions addressed, 
as we now explain. 

The main difference in the model we and [9, 10] consider is in our assump- 
tion (A) in contrast to their non-demolition assumption on the evolution. 
The latter assumes that there is a preferred basis of system states, called 
pointer states, which is preserved by the interacting system-probe dynamics. 
The pointer states remain unchanged under the successive measurements and 
evolution. As a consequence of this assumption, the system-probe interaction 
operator commutes with the system Hamiltonian, so there is no energy trans- 
fer between the system and the probes. Assuming this manifold of invariant 
states, it is shown in [9] that any initial state of the system converges, under 
the repeated measurement evolution, to one of the pointer states. Moreover, 
the measurement outcome determines which pointer state is chosen. The 
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results are derived for a homogeneous model (same setting for each inter- 
action cycle) and under a 'non- degeneracy condition'. They are generalized 
in [10] to non-homogeneous settings and without the non-degeneracy con- 
dition (under which the final dynamics becomes more complicated). Our 
assumption (A) is in some sense exactly the 'opposite' of a non-demolition 
assumption. Namely, it forces the dynamics (without measurement) to have 
a single stationary state (as opposed to an entire basis of stationary states). 
Both assumptions are reasonable, but they discribe different physical situa- 
tions. As pointed out in [9, 10], the non-demolition setting is realized in some 
experiments in quantum optics. Our setting describes scattering processes 
where energy is exchanged. A typical example is that of a 'one atom maser', 
where atoms (probes) interact with modes of the electromagnetic field in a 
cavity (system) by exciting the field modes, leading to subsequent photon 
emission [27]. These processes necessitate energy exchange. A famous model 
describing this situation is the Jaynes-Cummings model, in wich energy is 
not conserved, and which we discuss in some detail in Section 3. 

According to the different physical setups, the questions addressed here 
and in [9, 10] are different. Their goal is to gain knowledge on the system by 
measuring the outcoming probes. This can be done since both the system 
and the measurements converge for large times. (The convergence of the 
measurement outcomes is not mentioned as one of the main results in [9, 10], 
but it is established in the text of the article, see e.g. the ir^i) in [9].) 
As our results show, in the energy-exchanging situation, the measurement 
process typically does not converge. We then focus on the question how, 
on average, the input (incoming probes) is related to the output (outgoing 
probes) of the scattering process. Related to this is that incoming probes 
are taken to be in pure states in [9, 10] (they may be selected randomly from 
a set of pure states), while in our work, they may be mixed. As we explain 
in the intoduction above, the statistical uncertainty in the mixed incoming 
states contributes to the fluctuation of the measurement process. The inverse 
scattering problem, namely statements about the scatterer (system), we have 
only addressed in the exceptional situation where the measurement process 
does indeed converge (see Lemma 13). It could be interesting to study suit- 
ably averaged (and hence convergent) properties of the system state. This 
issue will have to be addressed elsewhere. 

On top of generalizing the results of [9] to a more general non-demolition 
setting (with random characteristics of the system, and e.g. possible feed- 
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back in the measurement process, and also considering the continuous-time 
limit), the authors in [10] prove the (almost sure) convergence of measure- 
ment outcome frequencies (sections 4.1 and 4.3 in [10]). This is analogous 
to the convergence statement in our Theorem 4. Their argument is based 
on the law of large numbers for Markov chain states. Our approach gives a 
concrete construction of the final outcome frequencies which is amenable to 
explicit calculations (analytic perturbation theory) and physical interpreta- 
tion (Theorem 4). 

There is an interesting similarity between repeated interaction and quan- 
tum random walk models. The literature on the latter is huge and it is not 
our goal to present it here, with the exception of [2], which has been pointed 
out to us by the referee of the present paper. In that reference, the authors 
consider the dynamics of an open system quantum random walker. The pair 
(p n ,X n ) defines a quantum trajectory process, where p n is the state of the 
walker and X n is its position (in Z d ). Under a suitable ergodicity assump- 
tion on the dymamics of p n (resembling our condition (A)), a central limit 
theorem is proven for the random process X n . In a sense, the X n is anal- 
ogous to our measurement outcome process, and it would be interesting to 
establish a central limit theorem in our case as well, giving a refinement of 
the law or large numbers established in Theorem 5. This question, together 
with an elucidation of the relation between repeated interaction systems and 
quantum random walks, presents a circle of interesting problems for future 
analysis. 

2 Quantum dynamical system setup 
2.1 Formalism 

The general formalism of quantum dynamical systems is presented in [5, 16] 
(see also [12]). Both S and V are described as quantum (W*) dynamical 
systems in standard form. The states of such a system are given by unit 
vectors in a Hilbert space "H, observables form a von Neumann algebra 9JT C 
B{l-L) and the dynamics is given by a group of * automorphisms a f on 9Jt, 
t £ i There is a distinguished (reference) vector ip e H which is cyclic 
and separating for 9Jt, and such that the dynamics is represented as a 1 (A) = 
e ltL Ae~ ltL , for A e 9JT, and where L is a selfadjoint operator on V. satisfying 
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Lip = 0. This operator is called the standard Liouville operator. 

Accordingly, the system S is determined by a Hilbert space Hs, a von 
Neumann algebra Wis, a cyclic and separating vector ip s and a dynamics 
a s = e ltLs ■ e~ ltLs , with L s ips = 0. A single probe is described by the 
same ingredients (with index S replaced by V). We assume throughout the 
paper that dim'Hs < oo, dim'H-p < oo. The Hilbert space of the chain of all 
probes is the tensor product He = ® n >\Hv, stabilized on the reference vector 
ip-p G Hp. The von Neumann algebra of observables of C is Wlc = ® n >\Wlv 
and its dynamics is a c = <8> n >ia:p. 

The full system is described by the Hilbert space 

U = Us®U c (2) 
on which acts the von Neumann algebra of observables 

m = Wl s ® Wl c . (3) 
The non-interacting Liouville operator is given by 

L = L s + J2 L n,r, (4) 

n>l 

where L n p is the operator acting trivially on all factors of H except on the 
n-th factor of He, on which it acts as L-p. The interaction between S and 
the n-th probe lasts for a duration of r > and is determined by an operator 
V n , acting trivially on all factors of "H except Hs and the n-th one in He, 
where it acts as a fixed selfadjoint operator 

V = V* G Wl s ® Wl v . (5) 

Let G "H be an initial state of the full system. Then the state of the 
system at time-step n is given by the vector 

*„ = */„••• Wi*o, (6) 

where 

U n = e iT ( Lo+v ^ (7) 
is the unitray generating the one-step time evolution at instant n. 
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We consider initial states of the form 

= ^5 ®n>i Bifo>, (8) 

where B is an operator in the commutant von Neumann algebra 

Wl' e = {A G B(U V ) : AX = XA \/X G Tic}, 

such that H-B^-pll = 1 and LpB = BL-p (normalized, invariant state). Since 
ip-p is cyclic for VJl'p separating for OJt-p, where 9Jtp the commutant of OJt-p) 
and since dim Hp is finite, every tp G H-p is exactly represented as ip = Bipp 
for a unique B G WX' V . 

Let J and A denote the the modular conjugation and the modular oper- 
ator associated to the cyclic and separating vector ips ® ipp for the von Neu- 
mann algebra Wls^Wlp. By the Tomita-Takesaki theorem [8], we know that 
A it (DJl s (g)Dnp)A- it = Wl s (g)DJlp for all t G R (A _i * is a group of unitaries), and 
that J(Wl s <g> J = 9Jt^ <g> £0^,. Consequently, JA^FA -1 * J G 371^ <g> 97^>, 
where V is the interaction operator defined in (5). In the finite-dimensional 
case as considered here, an easy analyticity argument shows that the last 
relation stays valid for any t G C. In particular, 

JA 1/2 VA- 1/2 J eWl' s ® m'p. (9) 

It will be convenient to represent the joint dynamics of the system S 
and the probe V interacting at the given moment with S by the following 
operator acting on Hs <8> Tip (see Subsection 2.3), 

K = L s + L v + XV - XJA 1/2 VA- 1/2 J. (10) 
Here, A 6 R is a coupling constant. Due to property (9) we have 

e itK Ae~ itK = e H i L s+Lv+w} Ae -it[L s +L P +xv] ^-q 

for all A G 9Rs <S> VJlp and all t G R, as is not hard to see for instance by 
using the Trotter product formula. The term —\JA 1 l 2 VA~ l / 2 J in (10) is 
introduced in order to have the property 

Kfc ® = 0. (12) 

The latter relation follows from (Ls + Lp)ips <8> tyv — and the fact that 
JA l / 2 V A" 1 / 2 Jip s ® i>v — Vips ® V'P (which in turn is implied by A -1 / 2 J = 
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JA 1 / 2 and JA^Atps ® ip P = A*ip s <g> ip P for all A e Tl s ® Wl v ). The 
operator K has been used in [25] for the study of non-equilibrium open 
quantum systems and in the setting of repeated interaction open systems in 
[12, 13, 14, 15]. 

2.2 Multitime measurement process 

In this subsection, we describe the process of multitime measurement of the 
outcoming probes. We refer to [29, 17, 24, 19] for a detailed introduction to 
quantum measurement theory. Let M e DJl-p be a selfadjoint "measurement" 
operator on H-p with spectrum spec(M) = {mi, . . . , m M }, where 1 < fj, < 
dim "Hp (distinct eigenvalues). Let S be any subset of spec(M) and denote 
by E s the spectral projection of M associated to S. 

Suppose that the entire system is in a state ^ EH initially, and that the 
following experiment is performed: the system evolves according to U\ and 
then a measurement of the observable M is made on the outcoming probe, 
yielding a value in S 1 C spec(M), and then the system evolves according to 
U 2 and after this evolution a measurement of M is made on the outcoming 
probe and yields a result lying in S 2 C spec(M), and this procedure is re- 
peated n times. According to the principles of quantum mechanics, the prob- 
ability for obtaining the multitime measurement result lying in Si, . . . , S n is 
given by 

P(Si, ...,S n ) = \\E Sn U n ■ ■ ■ E S2 U 2 E Sl U^ \\ 2 ■ (13) 

Furthermore, if this probability is nonzero (so that the outcome of the specific 
experiment is actually realizable), then the state of the system immediately 
after the n-th measurement is given by the normalized vector 

^ = E Sn U n ---E S2 U 2 E Sl U^ 
a/ P(Si, . . . , S n ) 

We have -E sp ec(M) = 1, which corresponds to the situation where at the given 
time step no measurement is performed. 

The stochastic process associated to the measurements is constructed as 
follows. Let 

Q = S N = {u = (ui, u 2 , . . .) : UJj e spec(M)} 
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and let J 7 be the a-algebra of subsets of Vt generated by all cylinder sets of 
the form 

{wGfl : tdi G Si, . . . , o> n G S n , fiGN, Sj C spec(M)}. 

On (fl, J 7 ) we define the random variables X n : Q — > spec(M) by X n (uu) = 
u n , for n = 1,2,... The random variable X n represents the outcome of 
the measurement at time-step n. The finite-dimensional distribution of the 
process {X n } n >! is given by 

P(Xi G Si, ... , X n G S n ) = P(Si, S n ), (15) 

for any n G N, any subsets Si, . . . , S n of spec(M) and where the right hand 
side is defined in (13). P extends uniquely to a probability measure on (f2, J 7 ) 
by the Kolmogorov extension theorem. 



2.3 Representation of joint probabilities 

The main result of this subsection is formula (22). The expectation of an ob- 
servable A G DJts of S at time step n is given by (A) n = (fy n , A$ n ), where 
defined in (14), is the normalization of the vector 4> n = E Sn U n ■ ■ -EsJJi^q. 
Note that only the n first factors in the product of the probe Hilbert space 
%c ( see before (2)) play a nontrivial role in the definition of \l/ n . We have 



*n,A* n ) (16) 

= (^o, e iT ^E Sl ■ ■ ■ e irZ -E Sn A E Sn ^ T ~ L - ■ ■ ■ Eg^^o) 

= (* , lU+Ye^E Sl ■ ■■e iTL -E Sn [U-}*AU-E Sn e- iTL - ■ ■ ■ ^e"^ 1 



where 

L^Ls + L^ + Vj (17) 

acts non-trivially only on the Hilbert space of S and the j-th probe Hilbert 
space, and where we introduced the unitaries 



U n =exp 



-irj^(n- j)L jtV 



and = ex P 



3=1 



(18) 
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Since A commutes with we have [U~]*AU~ = A. The operator Cj = 
e -"'0'- 1 ) i 3,p5 e 11 "(i- 1 ) i i,p belongs to the commutant Wl' p , since B does, and 
since the dynamics generated by L^-p leaves DJi' v invariant. Recalling that 
= (I5 ® B <g) B ■ ■ ■ ^rcf, with \l/ rc f = ip s ® i>v ® ' " j we obtain from 
(16) 

(y n ,AV n ^ 

= (d • • • C n * ref , e iril £ 5l • • • e iTi "£ s „A£ s „e- iTt " • • • Es.e'^C, • • • C n V rei ) 
= (* ref , CTd • • • C;C n e Ml i? Sl • • • j TKn E Sn A* ref ) 

= (* ref , [PxCTCie^^Pi] • • • [P n C* n Cj TK -E Sn P n \A^^) , (19) 

where Pj is the projection acting trivially on all factors of H except on the 
j-th "H-p, where it acts as the rank-one orthogonal projection onto ip-p. We 
define the operator Tj = PjC*Cje 1TKj EsjPj, which, under the hypothesis 
that e isL vBe~ isLv = B, becomes Tj = PjB*Be irK i E Sj Pj. We identify Tj as 
an operator on the Hilbert space Us, and as such, write 

T s = PB*Be irK E s P, (20) 

where P the orthogonal projection onto ip-p £g> 7-Lc, S C spec(M) determines 
the measurement performed at the given time-step, and where K is given in 
(10). We write simply T for T spec (M)- Remark that since Esj G VR-p and tp-p 
is separating for Tip, we have E Sj Pj 7^ for all j. With this definition, we 
arrive at 

(* n , At> n ) = ty s , T Sl T S2 ■ ■ ■ T Sn A^ s ) . (21) 
In particular, the probability (14) can be expressed as 

P(S 1 ,...,S n ) = {il>s,T Sl ---T Sn rl> s ). (22) 

2.4 Analysis of joint probabilities 

Lemma 6 The spectrum ofT$, (20), lies in the closed unit disk centered at 
the origin of the complex plane. For S = spec(M), i.e., E s = H, we have in 
addition Tip s = ip s . 

We consider the probability P(X n e S eventually), for S C spec(M). 
This quantity can be expressed using the Riesz spectral projections II and 
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lis of the operators T and Ts associated to the eigenvalue 1. They are defined 
by 

u s = ^-^(z-Tsy'dz, n = n spec(M) , (23) 

where the integral is over a simple closed contour in the complex plane en- 
circling no spectrum of T# except the point 1. If 1 is not an eigenvalue then 
II5 = 0. For the next result, we recall the following definition, 

{X n G S eventually } = {u\ there exists a k s.t. X n (oj) G S for all n > k}. 

Lemma 7 We have P(X n G S eventually) = (ips, n II5 ips). 

Proof. The set {X n G S eventually} is the increasing union of {X n G 
S Mn> k}, so P(X n G S eventually) = hm k ^ ao P(X n G S Mn > k). Next, 
{X n G S Vn > k} is the intersection of the decreasing sequence {X n <E S, n = 
k, . . . , k + /}, so 

PpT n G 5 eventually) = lim lim P(X n e S, n = k, . . . ,k + I). (24) 

fc— >oo I— >oo 

We have P(X n e S, n — k,...,k + l) — (ip s , T^T 1 / 1 ^) ■ Since for each 
k fixed, the limit of {^s-,T k ~ 1 T l g' l ^s) as I — > 00 exists (it is the probability 
P{X n e S\/n> k)), we have 

1 L 

lim (ij s , T k ~ l T l s +l ^ s ) = hm - ^ <^ 5 , T*" 1 !* 1 " 1 ^) = T*"^^) , 

00 00 ^ 

where II5 is the ergodic projection of Ts associated to the eigenvalue 1. 
(lis — if 1 is not an eigenvalue of Ts.) Arguing in the same way for the 

limit k — > 00, we obtain P{X n G S eventually) = ^SjIHIs^s^. Invoking 

Lemma 20 we replace the ergodic projections with the Riesz projections. ■ 

The probability is given entirely by information on the spectrum of T and 
T s at the point 1 (Riesz projections). The result holds even if T or T s have 
spectrum on the unit cirlce other than possibly at 1. 

Finer information about the asymptotic dynamics depends on all the spec- 
trum of T on the unit circle. We make the following ergodicity assumption 
(compare with Lemma 6). 
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Condition A. The point z = 1 is the unique eigenvalue of T with \z\ = 1, 
and this eigenvalue is simple (with eigenvector ips)- We define the gap by 

7 = l-sup-{>| : z G spec(T), z ^ 1}. (25) 

For A = the operator T is e 1TLs and has spectrum on the unit circle, 
with degenerate eigenvalue 1. Assumption A is verified in practice typically 
by perturbation theory (A small, nonzero). It is sometimes called a "Fermi 
golden rule condition". In this setting, condition A implies 7 > for small 
nonzero A. Condition A implies the dynamical behaviour of the assumption 
(A) stated before Theorem 1. 

The random variables X n are not independent but their correlations de- 
cay. 

Theorem 8 (Decay of correlations) Suppose Condition A holds. For any 
e > there is a constant C e such that for \<k<l<m<n< 00, 
A G o-(X k , . . . Xi) and B G cr(X m , . . . , X n ), we have 

\P(A f]B)- P(A)P(B)\ < C e P(A) e -(— Olind^)-^ ( 26 ) 

We give a proof of Theorem 8 in Section 4.1. For fixed e < ln( jz^), define 
the function C : N ->■ R+ by 

C(d) = C e e' dM ^ ) - e] . (27) 

Theorem 8 implies that the random variables X n have decaying correlations 
in the following sense: for all 1 < k < I < m < n, all A G cr(X k , . . . X{) and 
all B G cr(X m , . . . , X n ), we have 

\P(A(lB)-P(A)P(B)\<C(m-l), (28) 

where the function C is independent of A, B, k, I, m, n and satisfies C(d) — > 
as d — > 00. 

Let Xk, k = 1,2, ... be a sequence of random variables. We denote by 
a(X n , X n+ i, . . .) the sigma algebra generated by {Xk}k>n- The tail sigma 
algebra of the process {X n } n >i is defined by T = P\ n >io(X n , X n+ i, . . .). 

The following result is a generalization of the Kolmogorov zero-one law, 
valid for a process {X n } where the random variables are not independent, 
but have decaying correlations. 
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Theorem 9 (Extended Kolmogorov zero-one law) Let X n be a sequen- 
ce of random variables with decaying correlations, as in (28). Then we have 
P(A) = or P(A) = 1 for any tail event AeT. 

A proof of this result can be obtained by extending proofs of the Kol- 
mogorov zero-one law for independent variables, see [1]. Under Condition A, 
we write the rank-one Riesz projection of T associated to z = 1 as 

n = fofe)<^l, (29) 

where T^ s = ij) S , T*tp* s = tp* s and (ips^s) = 1, H^sll = 1. (T* is the adjoint 
ofT.) 

Corollary 10 Assume that Condition A holds and let S C spec(M). Then 

P(X n e S eventually) = (ip* s ,U s ips) e {0, 1}. 

Remark. Both and IL5 depend on A (ips does not). If condition 
A holds for A G /\{0} for some neighbourhood I C M. of zero, then for A 
sufficiently small we have P(X n e S eventually) = lim A ^ (^sW^sW^s)- 
This follows from the facts that the mapping A 1— > (ipg,Usips) is continuous 
in a deleted neighbourhood of A = and that the image is discrete. We point 
out that the map A i-> (ipg,Usips) is actually holomorphic in a punctured 
neighbourhood of the origin, and stays bounded there (the only possible 
image points being or 1, even for complex A, by the identity principle). 
Hence zero is a removable singularity of this map. As objects on their own, 
ipg and lis are not holomorphic at the origin in general (eigenvalue splitting), 
but their combination as in the inner product is. 

Theorem 3 gives a criterion for P(X n e S eventually) = 0. The next 
result is a characterization of when this probability is one. 

Lemma 11 Assume that Condition A holds. If (ips^s^s) 7^ for some 
small enough A then P(X n e S eventually) = 1 for sufficiently small values 
ofX. 

Proof of Lemma 11. The inequality P(X n e S eventually) > P(X n e 
S Vn > 1) is the same as (^5, IT Hs^s) > {^s, H-s^s) (see also Appendix A). 
Result (b) follows from the fact that the left side can only take the values 
zero or one, independently of A for A sufficiently small (guaranteeing that 
(ips, n Hsips) is continuous). ■ 
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Proposition 12 Let Aj G cr(Xj), j > 1. We /iawe 

sup |P(A n , . . . , A n+fc ) - P(A„) • • • P(A n+k ) | < C fc 

n>l 

for any k > 1, and for some constant C k . 

Proof of Proposition 12. It suffices to show that 

p(x n e s n , . . . , x n+fc g sVt+fc) — P(x n g s n ) ■ ■ ■ p{x n+k g s n+ k) = o(||v||), 

uniformly in n, and where Sj C spec(M). Since T 5 = PB*Be [rK E s P = 
e irLs u in (E s ) + 0(||V||), the joint probability on the left side is 

(■0s, T n ~ l T Sn ■ --Ts^s) = ujUEsJ ■ ■■"in(E Sn+k ) + 0(\\V\\). 

Similarly, P(Xj G Sj) = u in (E Sj ) + 0(||V||), and so the result follows. ■ 

In the next result we investigate the state of S knowing that X n G S 
eventually. Let u n be the state of S at time step n (obtained by reducing 
the state of the entire system to S). 

Lemma 13 (Asymptotic state of S) Suppose that the measurement out- 
comes 

X\ G Si, ... , X n -i G S n -i, X k ES,k>n (30) 

are observed for some n > 1, and that 1 is a simple eigenvalue of T s with 
Riesz projection lis — We have for any observable A of S 

{r,Ms) 



lim co n (A) = uJooiA) = 



This is an inverse scattering result: knowing that the scattered particles 
are measured to lie in S we can deduce the state of the scattering object S. 
The final state does not depend on the initial outcomes X n for n < k, any k. 
However we show in the proof of Lemma 13 that if the eigenvalue 1 of Ts is 
not simple, then the system converges to a final state which depends on the 
whole measurement path X±, X 2 , . . . 

Proof of Lemma 13. The asymptotic state of the system is 

, A s y (^ S ,T 1 ---T n . 1 T s Ajj s ) 

iOooiA) = hm —. : 

{ips,T 1 ---T n - 1 T l s ip s ) 
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see (21). Now (ips,Ti ■ ■ -Tn-iTgips) converges to the nonzero probability of 
observing (30). Therefore, by (the proof of) Lemma 20 

lim (fa, 7\ • • • T^T&s) = U s , T x • • • T n _ 1 n^ (S ) = ($ s , T\ • • • T n _^ s ) ■ 

By simplicity of the eigenvalue 1 of T s , 



^oo (-4) 



(^ 5 ,T 1 ---T n _ 1 ^)(^,^) ' 



If 1 is not a simple eigenvalue of T s , so that Tl s — Y7j=i iV'j)^!) the final 
state is _ 

_ E; = i (jgg, T, ■ ■ ■ T n _iipj) (g, A^) 

The final state then depends on the whole measurement history. ■ 



3 The Jaynes-Cummings model 

We consider a simple system where both the scatterer and the probes have 
only two degrees of freedom participating in the scattering process. The pure 
state space of S and V is C 2 , and the Hamiltonians are given by the Pauli 
<t 2 operator, 



He = H% 



1 
-1 



(31) 



The interaction between S and V is determined by the operator 

XV = A (a% <S> a v + a s ® a* v ) , (32) 
with coupling constant A G 1, and where 





1 



a = 



1 




(33) 



are the annihilation and creation operators. In the usual Jaynes-Cummings 
model (used e.g. in quantum optics), the system S has many levels (harmonic 
oscillator), see e.g. [29] and also [18]. Our model is a truncation, but it still 
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describes energy exchange between S and V . In what follows, we can treat 
all values of A, not necessarily small ones only. This is so since the model 
is essentially exactly solvable. The total Hamiltonian H = H s + H v + XV 
describes exchange of energy between S and V, while the total number of 
excitations, N = a^a^ + a^a-p, is conserved (commutes with H). This allows 
for a treatment of the system separately in the invariant sectors N — 0, 1, 2. 

For an arbitrary probe observable X G £>(C 2 ) we write Xij = ((pi,Xipj), 
where ipi, ip 2 are the orthonormal eigenvectors of H-p (with H-pifi = ipi). 
Incoming states are determined by p G [0, 1] via 

u- m (X)=pX 11 + (l-p)X 22 , (34) 

where X G B(C 2 ) is an arbitrary probe observable. 

We will use the notation and definitions of Section 2 in what follows. In 
particular, the single step operator Ts is defined in (20). For the following 
explicit formula, we take the reference state ^5 to be the trace state. 

Theorem 14 (Explicit reduced dynamics operator) Set ipij = ifi® ifj 

and let X be any operator ofV. In the basis {<fu, </?i2, <£>2i, ^22} we have 

PB*Be irK XP = u in (X) e [rLs + (35) 



(1 - p)X 22 a (1 - p)X 21 b 

-pX 12 e 2iT ism(\T) e 2ir (cos(Ar) - l)w m (X) 

pX 21 e- 2iT ism(\T) 

-pX 22 a -pX 2 ib 

-{l-p)X l2 b -(l-p)Xna 

(1 -p)Xi 2 e 2iT isin(Ar) 

e- 2iT (cos(Ar) - l)uj- m {X) -(1 - p)X 21 e- 2ir isin(Ar) 

pX\ 2 b pXna 



where a = — sin 2 (At) ; b = — isin(Ar) cos(At). 

We point out that the vector \p, 0, 0, 1 — pf is an eigenvector of the adjoint 
of (35) with eigenvalue Wj n (X). 

Proof. The proof is a rather lengthy but straightforward, explicit calcula- 
tion. Since H s and H v commute with I = V-V, where V = JA 1/2 VA~ 1/2 J 
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(see (10)), it suffices to calculate 



PB*Be iTl XP = J2 ( k ) {-l) n ~ k PB*BV k {V') n - k XP. 

n=0 ' k=0 ^ ' 

Here it is understood that all operators are considered in the "doubled" 
(GNS) Hilbert space, e.g., 

V = a* s <g> l s <g) a P <g> l v + a s ® I5 <E> Op <E> V 

Powers of V and V' can be calculated explicitly. For instance, for k > 2 even, 
we have V k = n k/2 ® 1 <g) (1 - n) fc/2 <g> 1 + (1 - n) fc/2 (g) 1 <g) n fe/2 <8> 1, where 
n = a*a. One obtains similar expressions for k odd, and for {V') 1 . Using 
these expression in the above series, one gets the result of Theorem 14. ■ 

Resonant and non-resonant system. If Ar is a multiple of 7r then 
(35) reduces to PB*Be lrK XP = u in (X) diag(l, ±1, ±1, 1) with plus and mi- 
nus signs if the multiple is even and odd, respectively. Then, by using the 
expression P{X± e Si,...,X n e S n ) given in (22), it is readily seen that 
the random variables Xj are independent, and P(Xj e S) — u- m (E s ). When 
Ar G 7rZ we call the system resonant [18], otherwise we call it non-resonant. 
One can understand the resonant regime as follows: consider the dynamics 
generated on S and a single probe V by the Hamiltonian H = H s + H-p + W . 
The probability of transition from the initial state <pf®<p±i where the S is in 
the ground state and V in the excited state, to the opposite state ff®f~2, at 
time t, is given by P t = \((pf <g> (p%,e itH (p$ <g> ifii)\ 2 = sin 2 (At). For At G 7rZ 
this probability vanishes. If the interaction time r in the repeated interaction 
system is a multiple of tt/X, then interaction effects are suppressed. It is not 
hard to see that in this case, the system does not feel the interaction with the 
probes in the sense that ui n (A) = uj (a^(A)) for all n > 1, where a^(A) is the 
reduced dynamics of S alone. We focus now on the non-resonant situation. 

Asymptotics. We show that 

For u>i n = |t)(t|; the process X n converges to |t)(tl (meaning thatP(X n = 
1 1) eventually) = 1) and the incoming state of the probes is copied onto the 
scatterer, in the limit of large times. 

This result is non-perturbative and holds for all A G R. We point out that 
this copying mechanism has been described before ("homogenization"), see 
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[11] and references therein. (In this respect, our analysis is more complete, 
as it describes the entire system of scatterer and probes.) Let oj- in = |t)(t|> 
corresponding to p — 1. The matrix in (35) becomes lower triangular, and 

spec(PB*Be irK XP) = X n {1, e 2ir cos(Ar), e~ 2ir cos(Ar), cos 2 (Ar)}. 

Condition A is always satisfied in the non-resonant case. ips and ipu are the 
eigenvectors associated to the eigenvalue 1 of PB*Be lrK P and its adjoint, 
respectively. So the Riesz projection associated with the eigenvalue 1 is 
LT = v / 2|'?/'5)(v 9 ii| ( see (23)). We can now readily obtain the following results, 
for an arbitrary measurement operator M. 

The asymptotic mean is given by /i^ = u in (M) (expression (61)). The 
frequencies are f m = w in (i? m ) (expression (52)). This suggests that the cavity 
becomes 'transparent' for large times (no effect on incoming probes). And 
indeed, the Riesz projection II5 associated to the value 1 of Ts (with S C 
spec(M)) vanishes unless (E s )n = 1. Thus we have P(X n e S eventually) = 
if (E s )n 1 (see Lemma 7). If (E s )u = 1 then E s is the projection 
and one finds easily that LT5 = \(pu)((p n \ and (using Lemma 7), 
that P(X n = I ^) eventually) = 1. In the situation where the measurement 
outcomes stabilize, X n — > we can determine the asymptotic state of the 
scatterer S from Lemma 13. It is given by Uoo(A) = u in (A), A e B(C 2 ). The 
above 'asymptotic transparency' statement follows. 

Large deviations for the mean. The logarithmic moment generating 
function [23] is defined by 

A(a) = lim ilogE[e naZ "], (36) 

for a 6 K si. the limit exists as an extended real number. Using expression 
(66) and Theorem 14 (with p = 1), we find that A(a) = logw in (e aM ), for 
a e R. The Legendre transformation of A(a), 

A*(x) = sup ax - A(a), (37) 

x G R, is called the rate function. Its usefulness in the present context is 
due to the Gartner- Ellis theorem [23], which asserts that for any closed set 
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FcM and any open set G C 1, we have 

lim sup - log P (X n e F) < - inf A* (x) 

n^oo U x<=F 

lim inf — log P (X n e G) > - inf A*(x). (38) 

Here, T denotes the set of 'exposed points of A*'. A stronger "full large 
deviation principle" holds if an additional condition, called steepness of A, is 
verified. (See [23], Theorem 2.3.6 (c).) However, this condition does not hold 
in the present example. (Indeed, |A'(a)| -/^ oo as a — > ±oo, as is apparent 
from the explicit form A(a) = logu;i n (e aM ).) 

Proposition 15 Suppose thatVar(M) = u- m (M 2 ) — u- m (M) 2 , the variance of 
M in the state uj- m , does not vanish. Then A* is holomorphic at x = Ui Q (M), 
and 

Proof. Note that A is twice differentiable, and the second derivative 
w.r.t. a of the argument of the supermum in (37) is less than or equal to 
zero. Therefore, for fixed x, the supremum is taken at a G K. satisfying 

For a = we have x = u in (M). If A"(0) = Var(M) ^ 0, then equation (39) 
has an implicit solution a = a(x), locally around x = u in (M). Since A' (a) 
is holomorphic at a — 0, the implicit solution is holomorphic at x — w in (M) 
(see e.g. [32], p. 163, equation (12.4)). The Taylor expansion of (39) is 

x = u in {M) + aV&i{M)+ca 2 + 0{a 3 ), (40) 

where c = |{w in (M 3 ) - 3w in (M 2 )w in (M) + 2w in (M) 3 }. We solve equation 
(40) implicitly for a = a(x), which is the point where the supremum in (37) 
is taken. The explicit formula for the supremum given in Proposition 15 
follows. ■ 

Example: Measuring the outgoing spin angle. Since u- m is the state 'spin 
up', we have u in (M) = M u and Var(M) = |M 12 | 2 . Imagine an experiment 
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where we measure the angle of the spins as they exit the scattering process. 
Let 9 G [0, 7r/2) and G [0, 2n) be the angles measuring the altitude (0 = 
is spin up) and azimuth (0 = is the plane orthogonal to the axis of the 
cavity). The measurement operator "spin in direction (9, 0)" is given by 



M 



cos 9 e ^sm9 
e 1 * sin 9 — cos 9 



see e.g. [20], Chapitre IV, (A-19). The eigenvectors of M associated to the 
eigenvalues ±1 of M are 

X+ = e^' 2 cos(#/2)v?i + e^ /2 sin(#/2)v? 2 
X- = -e -i * /2 sin(0/2)^i + e^ /2 cos(0/2)y9 2 . 

The eigenprojection E + measures the spin in the positive direction (9, 0). By 
using Lemma 7 is easy to see that 



P (X n is in direction (9, 0) eventually) 



1 if = 
if 0^0. 

This is another manifestation of the asymptotic transparency of the cavity. 
We obtain from theorem 5 (with /i^ = cos 9) that for any e > 0, 

lim P(\X n - cos0| > e) = 0. 

n— ¥oo 

The speed of convergence can be estimated using (38) and Proposition 15. It 
is easy to see that the logarithmic generating function and the rate function 
associated to the shifted random variable X n — cos 9 are given by A S hift(o;) = 
A (a) — a cos 9 and A* hift (x) = A*(x + cos 9), respectively. Next, we note that 
all points in the vicinity of zero belong to the set J^hift, the set of exposed 
points of A* hift . Indeed, if x = A' shift (a) for some ael, then x G J^hift ([23], 
Lemma 2.3.9). But x = = Ag hift (0), and A' shift is invertible around zero (as 
A" h;ft (0) 7^ 0). This shows that J^hift contains a neighbourhood of the origin. 

Take < e < e' « 1, set G — (— e', — e)U(e, e'), and let F be the closure of 
G. Then mi xeF A* shiit (x) = mi xeGnTshm A* shiit (x) = + 0((e') 4 )- (We use 

Proposition 15.) Combining this with the two bounds (38) (for the shifted 
random variable), we obtain 

r 2 



P(e < \X n - cos0| < e') ~ exp 



-"{d?i + " e '» 4 »} 



n — > oo, 



which is a large deviation statement for the average X n . 
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4 Proofs 



4.1 Proof of Theorems 1 and 8 

Theorem 8 is a stronger version of Theorem 1, so it suffices to prove the 
former. 

Let A G cr(X k , . . . , X/) and B G cr(X m , . . . , X n ). The range £ of the X n 
is finite, so a(X k , . . . ,Xi) consists of the collection of all sets of the form 
{u : (X k (co), . . .,Xi(u)) G H}, where H C £'- fe+1 ([8], Thm. 5.1). Therefore 
we have 



A = (jxt\{8?})n...nxr\{s^)=:\jA j 

3=1 3=1 

b = U««)n---nV({ 8 W})=:UB il 

1=1 1=1 
where G £, and ^ n Aj = 0, P^ n P, = for i ^ j. Thus, 

p(A) = £p(a,-), P(s) = 5]p(s i ), P(AnB) = J £p(A j nB i ). 

j i i,3 

Setting Tm^ := T , („>, , we have 

P(A nB) = ^<#s, T k ~ l T^ ■ ■ ■ T^T m ~ l ~ x T^ ■ ■ ■ rW* s >. (41) 



We now approximate T m ~ l ~ x by its value for large m — I. To do so, let Ps 
denote the Riesz spectral rank-one projection onto the eigenvalue one of T, 
and let P s = l s — P s . We have P s = IV'sKV'sIj where ip* s G %s satisfies 
T*i>* s = ipg and ^5) = 1. The operator T has the spectral representation 
[26], I§5 



d 



T= mm +J2i Z r P r + Dr}, (42) 



= 1 



where P r is the Riesz projection associated to the eigenvalue Zj and D r is 
the associated eigen-nilpotent. We have D v r r = (u r is the index of z r ). 
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Consequently, 

d v r — l 



T k = Mini + E E 4~ q PrD q r =: |V5>(^I + Rk- (43) 

r=l g =0 



Note that it suffices to consider the nonzero z r in (43). Using the bound 

yv—l /k\ 
'9=0 \q) 



Ylq=o — V ^ V 1 — k u > an d the f &c ^ that f° r an y e > there is a constant 



C e s.t. A;^ < C t e ek for all fc > 1, we obtain 

\\R k \\ < C e e ek max \z r \ k < C e e ek {l - 7 ) fc , (44) 

l<r<d 

where we invoke Assumption A, \z r \ < 1 — 7. We now replace T™" 1 " 1 in (41) 
using (43), 

p(Ans) = p(A)5>;,rW...r£-ty fl ) (45) 

+ E(^5, T k ~ 1 T^ ■ ■ ■ T^R^tV ■ ■ ■ T^s) (46) 
The sum in (45) equals 

E^ T ™--- T ^) = E^^d^x^i)^---^) 

j i 

= E<V>5, (^ m_1 - • • • rw^> 

i 

= P(B)-J2^s,Rm-iT^---T^ s ). (47) 

i 

Lemma 16 There is a constant C s.t. \\ ^2 { Tm ■ ■ - T^\\ < C ; independently 
of the range of values of i, m, n and of the defining the 

We give a proof below. Lemma 16 together with (45), (46) and (47) gives 
\P(A n B) - P(A)P(B)\ < CP(A)\\R m ^\\ 



+ 



E(^, T k ~'T^ ■ ■ ■ T^R^T^ ■ ■ ■ T^ s ) 



h3 



(48) 
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The remainder i? m _/_i is given by the sum in (43) (with k — m — I — 1). We 
expand 

P r D« r = J2^ s ')\Xs)(xA (49) 



in an orthonormal basis {Xs}, where ou(s,s') G C are matrix elements (also 
depending on r, q). The modulus of the sum in (48) is bounded above by 



d u r —l 
r=l q=0 



m — I — 1 



U}(S, s 



X 



(50) 



Lemma 17 Lei A s G 971s fre £/ie unique operator s.t. Xs = A s ijjs. We /iai>e 



< Ps||P(A). 



We give a proof of Lemma 17 below. Using the result of the Lemma in (50) 
we obtain, for any e > 0, the upper bound P(A)C e e € ( m ~ l \l — 7)™"' for the 
sum in (48). Thus for all e > there is a C e s.t. 

\P(AnB)-P(A)P(B)\ < C e P(A){(l- 7 ) m e em + (l-7) m - / e e(m - /) }. 

Since m4(l- 7) m e em i s decreasing (for e small enough), we get (26). This 
concludes the proof of Theorem 8 modulo proofs of Lemmata 16 and 17. 

Proof of Lemma 16. Any vector ip G Us is of the form ip = A^/ s , for 
some A G 97t,s. We have 



rrt n <-> 



sup 

15*511 = 1 



sup I (BV S <g> 5^p„ 

|B* S ||=1 



irL n 



ArL n 



QA x 



xe 



-irL r . 



-irL n 



where Q := Z)j £^({5™ }) ' ' -^({s^}) is a selfadjoint projection, ||Q|| = 1. 
Thus || Ei T « • • -TPAVsW < \\A\\. The result now follows from the uniform 
boundedness principle. This proves Lemma 16. 
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Proof of Lemma 1 7. By using the definition T£ 



Cj) 



where E { k j) = E k ({s ( k J> }), see (20), we see that 



,0") 



P k B*Be WK E ( k 3) P k , 



where = £ 



GO 77i(i)„irL, 



£f 'e' 



(i)\ 



S^-p (/ probes). To 



arrive at the form of we replace the action of e lrK in the operators T k ^ 
(see just above) by the action of the Liouville operators e lrife -e~ 1Tifc , see (11), 
(19). Since a positive linear functional on VJls is bounded, with norm equal 
to its value for the observable 11 ([16] Prop. 2.3.11), we have 



k 



T^A^s) < \\As\\ (^, T k ~ L T k ! ■ ■ • 7™ life 



k-lrp(j) 



,00 1 



Note that the scalar product on the r.h.s. is non-negative, as it is a proba- 
bility. It now follows that 

J2 I^T^Tf • • -If Xs > < \\A S \\ ^s,T k - l T^ ■ --T^s), 



the latter sum being P(A). 



4.2 Proof of Theorem 3 

We have P(X n e 5 eventually) = (^5, niLj^s), where II and II5 are the 
Riesz projections of T and T5 respectively, associated to the point 1, see 
Lemma 7. We know that T$ = Ui n (Es) + 0(|| V||). Since Eg is an orthogonal 
projection and u in is a state, u in (E s ) 7^ 1 means u in (E s ) < 1. Thus, for small 
enough V, T s does not have an eigenvalue at the point 1, and consequently 
U s = 0. ■ 



4.3 Proof of Theorem 4 

The proof is in two parts. First we show that the expectation of f m has 
the indicated limit. Then we upgrade this convergence to almost everywhere 
convergence, using the decay of corrleations, Theorem 1. 

Step 1. We show convergence of the expectation of the frequency, 

F m : = E[/ m ] = lim -E [number of k e {1, . . . , n) s.t. X k = m\. (51) 

rt— s-oo Tl 
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Let m G spec(M) be fixed and define x by x( x ) — 1 if x — m an d x( x ) = 
otherwise. Then the expectation in (51) equals 



E 



rai,...,ro„ 



P(X 1 = m 1 ,...,X n = m n ) 



= ^2 S = mi '- ■ • = m ' • • • ' Xn = m ^ 

3=1 m k ,kytj 
n 

3=1 

Therefore, since \ YTj=l T J ' — > II as n — > oo (the ergodic projection which 
equals the Riesz projection of T associated to one, see Appendix A), we 
obtain 

F rn = (^ s ,UT m ^ s ). (52) 

Recall that T m = PB*Be irK E m P, where K = L + XI. To expand (52) in a 
power series in A, we use 

FT fTn-1 

e irK = e irL ^i\) n / d n ■ ■ • / dr n /(r 1 ) • • • J(r„), (53) 

n>0 J ° J ° 

where J(s) = e-' lsL °Ie isL ° . The series converges in operator norm, for all 
A G R. This gives T = e irLs + T', where 

T' = e iTLs V(iA)" dn • • • p 1 dr n PB*BI(r 1 ) ■ ■■I(r n )P (54) 



For zeCin the resolvent sets of e lrL5 and T, we have 

(T - z)- 1 = (e irLs - z)- 1 J2 {-T'(e irLs - z)- 1 )" . 



(55) 



n>0 



The last series converges for A small enough (T" is of order A). Using (55) 
and (54) in the expression II = ^[§(T — z)~ 1 dz, where the integral is on 
a small circle encircling the point one of the complex plane, but no other 
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spectrum of T, we obtain 
ITT 



= 2^if dz ^ TLS ~ 1 ~ iAgirLS / dT 1 PB*BI(T 1 )P(e irLs - z)' 1 

+X 2 e [rLs [ dri / 1 dT 2 PB*BI(r 1 )I(r 2 )P(e iTLs - 
Jo Jo 

- A 2 (V rLs jf dT 1 PB*BI(n)P(e irLs - z)- 1 ^ + 0(A 3 



xPB*B 



1 + iX dr 1 /(r 1 ) - A 2 / dn / dr 2 /(r 1 )/(r 2 ) + 0(A 3 ) 



P m P 



By taking averages in the state ifts, we obtain F m . The term of order A of 
F m is thus (use that L s ^>s = 0) 

1 



2m 



dz (ip s , 



e lrLs - z] 



- 1 e 1TLs PB*BE m Pij s ) = ou in (E m ). 



There are, at first sight, two terms of order A 1 . However, the one contain- 
ing two resolvents (e lrL<s — z)^ 1 vanishes, as by taking the average in the state 
ips, each resolvent reduces to (1 — z)^ 1 , and ^ § dz{\ — z)~ 2 = 0. The non- 
vanishing term of order A 1 in F m is iA (ips, PB*Be~ 1TlL °Ie 1TlL °E m Pilj s ^ dri. 
Writing I — V — J A l ' 2 V J A 1 ' 2 , for V E Wl s (g) DJl P , the integrand of the last 
integral has the form 

(^ s ,PB*B{V(r 1 ) - JA 1 / 2 \/(r 1 )JA 1 / 2 }P m P^ 5 ) 

= (fa ® fo, B*B{V(n) - JA^V^JA^jE^s ® i> v ) 
= (x/js ® ifr, B*B{V(n)E m - E m V(n)Ws ® i>v) 
= u s ®u- m {[V{T 1 ),E m }) . 

The form of f m given in Theorem 4 follows. We point out that in the above 
four formulas, we use the same notation for the operator V and its repre- 
sentation on the GNS space. This should not lead to any confusion, as the 
context makes it clear which operator we mean (and it has the advantage that 
we do not have to introduce further notation for the representation map). 
This finishes the first step of the proof, showing that 



lim — E [number of k e {1, . . . , n) s.t. X k = m] = u, (56) 

n— s-oo Ti 
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where v := Ui n (E m ) + \f' m + 0{\ 2 ) and f' m is defined in Theorem 4. 

Step 2. We upgrade (56) to almost sure convergence, by using a classical 
fourth moment method. Introduce the random variable 

Z n = {number of k e {1, . . . , n) s.t. X k = m}. 

We are going to show below that, for any e > 0, 

£ P(\Z n /n - E[Z n ]/n\ 2 > e) < oo. (57) 

n>l 

Then by the (first) Borel-Cantelli lemma, P(\Z n /n — E[Z n ]/n\ 2 > e i.o.) = 0, 
i.e., there is a set £ of measure one, s.t. for all w 6 E, there exists a A; with 

|Z n /n - E[Z n ]/n\ 2 < e, Vn > A;. (58) 

From (56) we know that E[Z n ]/n converges to z/, so (58) implies that Z n /n 
converges to v almost everywhere. It remains to prove the summability (57). 
By Chebyshev's inequality, 

P(\Z n /n - E[Z n ]/n\ 2 > e) < -^E[\Z n - E[Z n ]\ 4 }. (59) 
We get an upper bound on the r.h.s. Set 

n n 

Z n - E[Z n ] = J2ix(X t = m) - P(X t = m)} =: ^F,, 
i=i i=i 

so that E[Yi\ = 0. Here, x(Xi — m) — 1 if Xi — m and x(Xi = m) = 
otherwise. We have 

n 

E[\Z n -E[Z n ]\ 4 ] < MYi\\Yi\\Yk\\Yi\]. (60) 

i,j,k,l=l 

The idea is to control the sum by using that if the indices i,j,k,l are far 
apart from each other, then the expectation is small due to the decay of 
correlations. Thus only a few terms in the sum contribute to its value. Let 
A > be a given integer "length" scale. All combinations of values of the 
four indices k, I belong to exactly one of the following five cases: (1) all 
indices lie inside an interval of length A, (2) three indices lie within A, the 
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fourth does not, (3) two pairs of indices are separated by more than A, but 
within each pair, the indices are apart at most by A, (4) one pair lies within 
A, the other two indices are apart from each other and from the close pair 
by more than A, (5) all four indices are apart from each other by at least A. 
Let rij, j = 1, . . . , 5, be the number of terms in the sum that satisfy cases (1) 
to (5) above. We have n x < 4!nA 3 , n 2 < 4!n 2 A 2 , n 3 < 4!n 2 A 2 , n 4 < 4!n 3 A, 
< 4!n 4 . Next, due to the separation of indices and Theorem 1, each term 
of case (5) is of the form 

EN^-llnM <E[|F i |]E[|F J |]E[|r fe |]E[|r,|] + 0(e-T' A ) =0(e"T' A ), 

as E[|Yj|] = for all i. Similarly, each term of cases (4), (3) and (2) have 
the same upper bound. Each term of case (1) is bounded above by one. We 
conclude that 

E[|Z n - E[Z n ]| 4 ] < C{n 2 + n 3 + n A + n 5 ) e -^ A + n v 

Choose now A = n a , with < a < 2/3. Then using the above bounds on rij, 
we have that n~ 4 E[|Z n — E[Z n ]| 4 ] is summable over n > 1, i.e., by (59), the 
inequality (57) holds. This completes the proof. ■ 

4.4 Proof of Theorem 5 

We prove that lim n ^oo E[X„ — fi^] = 0, with 

^ = (fc, UPB*Be WK MP^ s ) . (61) 

By proceeding as in step 2 of the proof of Theorem 4 above, one upgrades 
this to almost everywhere convergence, X n — > fj,^ a.e. Since 

E[X n ] = - ^2 ( m i H h m n )P(X 1 =m u ...,X n = m n ) 

mi,...,m„ 

= -J2(ij s ,T^ 1 PB*Be iTK MP^ s ), (62) 

and ^ YTj=i T^ 1 — > II = LT as n — > oo, we obtain the law of large numbers 
with (61). To obtain the expansion in A for A small, we find the Taylor series 
of LT = n(A) and e lrK ( x \ as in (55) and (53). Note that this calculation is 
the same as that carried out in Section 4.3 to calculate the expansion of F m , 
under the replacement of E m with M. The result follows. ■ 
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A Ergodic and Riesz projections 



Let V be a power bounded operator on a Hilbert space TL, i.e., such that 
\\V n \\ < M for a constant M independent of n G N. Set 7 = {y? G H : 
V> = y?} and ft = Ran(l - V). 

Lemma 18 We have H = J 7 + ft (closure) with J 7 fl ft = {0}. Moreover, 
the projection onto J 7 in this decomposition is the ergodic projection, II = 

Proo/. Clearly i £* =1 \/> = y? for all iV and all (/? G J. Also, 
F £n=i ^ n (! - v )1> = TrVty ~ yN ^] -» as iV oo, for any ^ G U. 
Thus, if y? G J= n ft then y? = £ Eli ->■ as N ->■ oo. This shows 
J 7 H ft = {0}. Similarly one shows that J 7 n ft = {0}: let tp G J 7 n ft Then 
V? = linifc (/9fc, with G ft. We have 

TV N N 

n=l n=l n=l 

The first equality holds since y? G J 7 . Since V is power bounded, the norm 
of the second sum on the r.h.s. is bounded above by M\\ip — ipk\\ for some 
M independent of N, k. The first sum on the r.h.s. converges to zero as 
N — > oo, since <fk G ft. Thus, upon taking first k, then N large enough, we 
see that y? = 0. This shows that J 7 fl ft = {0}. 

The equality H = J 7 + ft is equivalent to "H = J 7 + ft. We have 

T 1 - = Ran(l - V*), 7Z ± = {y? G U : V*<p = y?} 

Let y? G (J 7 + ft) ± . Then y? G J 7-1 flft -1 . However, V* is power bounded and 
thus, as above, J 7 - 1 fl ft- 1 = {0}. This shows that J 7 + ft is dense in "H. ■ 

We have {IT = <^ J 7 = {0}} <^ {1 is not an eigenvalue of V}. Assume 
that there is a neighbourhood U of 1 in the complex plane which does not 
contain any spectrum of V, except possibly the point 1. Let IT be the Riesz 
spectral projection, II = ^ § T {z — V) -1 d;z, where r C U is a simple closed 
curve encircling 1. II acts as the identity on J 7 , as does II. Let y? = (1 — V)x G 
ft. Then Tl(p = (1 — V)n%. The operator (1 — V)U is the eigen-nilpotent 
associated to the eigenvalue 1 of V. If 1 is a semisimple eigenvalue of V, then 
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Hip = for all tp G TZ, and by continuity, UTZ = {0}, which coincides with 
the action of II on TZ. If the eigenvalue 1 is not semisimple (the nilpotent is 
nonzero), then the action of II on TZ does not coincide with that of II. This 
shows the following result. 

Lemma 19 II = II if and only if 1 is a semisimple eigenvalue ofV (or if 1 
is not an eigenvalue ofV, in which case II = II = 0). 

Lemma 20 Suppose that dim'H < oo. If 1 is an eigenvalue of V then it is 
semisimple. We thus have U — II. 

Proof. Suppose 1 is an eigenvalue of V. It suffices to show that z h-» 
(z — V)^ 1 has a simple pole at z — 1. We have (z — V)^ 1 = + (z — 

V)-\l - II), where V = (1 - U)V(1 - ft) f Ran(1 _n } , so we only need to 
show that 1 ^ spec(V). Suppose that 1 G spec(K) and take ip G Ran(l — II) 
satisfying Vip = tp, \\ip\\ = 1. Since UV = VU we have Vip = V<p, so 
{y n — l)ip = for all n — 1, 2, . . . By applying ^2^=i an< ^ taking — > oo, 
we obtain (II — l)tp = 0. This is in contradiction to tp G Ran(l — II) with 
\\tp\\ = 1. ■ 



B Logarithmic moment generating function 

It is possible to give general conditions ensuring that the logarihmic moment 
generating function (36) exists. However, these conditions will be rather 
abstract (see before (66) below). They amount to knowing that a certain 
operator (R x (a) given in (64) below) has a unique eigenvalue of largest mod- 
ulus, with corresponding eigenprojection satisfying a non-vanishing overlap 
condition. We show in this section how perturbation theory can be applied to 
analyze the spectrum of the operator in question. The formulas established 
here can be used in the analysis of concrete systems. (An example being the 
Jaynes-Cummings model of section 3.) 

To calculate the logarithmic moment generating function, (36), we write 

E [ e a(x 1+ ...+x„)] = e a(mi+ - +m " ) P(X 1 = m l ,...,X n = m n ) 

mi,...,m„ 

= u(e aM r(4>,R x (a) n ^), (63) 
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where 

R x (a) = J2 eamT {m}=PB*Be^ K A(a)P (64) 



m 



*<<•> = ^s-y (65) 

The scaling of A (a) by l/u(e aM ) gives the normalization Ro(a) = e lrLs . The 
existence of the limit (36) is guaranteed if R\(a) is diagonalizable and has a 
unique eigenvalue p+(A, a) of largest modulus, and the corresponding Riesz 
projection P+(A, a) satisfies (tp,P + (\,a)ip) ^ 0. In this case (36) and (63) 
give that 

A(a) = logw(e aM ) + logp + (A,a). (66) 

In concrete examples, it is usually not possible to explicitly evaluate p + 
(however, it is possible in the example considered in Section 3 !), so pertur- 
bation theory is in order. We analyze the analyticity properties of R\(a) in 
A and a. Let Pj be the orthogonal spectral projections of M. For a G C we 
have 



ly* e 2mjRea||p.^ 

\\A(a)P\ " ' 



E^ mja \\P^v\ 



21 



For a G K there are constants < c < C < oo, independent of a G M, s.t. 
c < || A(a)P|| < C. This cannot be extended to all complex a, since otherwise 
A(a)P would have to be constant in a by Liouville's theorem of complex 
analysis. However, if the imaginary part of a is small, then the weighted 
superposition of the HP^-pH 2 of the denominator is still bounded away from 
zero. The growth of the numerator and denomiator as Rea — > ±oo is the 
same. Thus there is an a > s.t. if |Imo;| < a 0) then d < ||A(a)P|| < C 
for some < d < C < oo. 

The operator R\(a) is holomorphic in (A, a) G C x {z : \lmz\ < a }, 
and has the expansion 

P A (a) = ^A"pW(a), (67) 

n>l 

with 

p(«) ( a ) = i n d TLs [ dsi ■ ■ ■ [ " 1 ds n PB*BI(s n ) ■ ■ ■ I( Sl )A(a)P, (68) 
Jo Jo 
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and where /(s) = e isL °ie- isL °. The bound ||i2 (n) (a)|| < ^ff^- \\B*B\\ \\A(a)P\ 
implies that 

^11 T 1 1 ft 

sup ||i2 (n) (a)|| < C"—^- (69) 

|Ima|<ao n - 

for some C" < oo. Thanks to this bound we can perform perturbation theory 
in A uniformly in a s.t. |Imai| < a - 

Proposition 21 There are constants C , A i; both independent of a e C wii/i 
|Ima| < a and of t > 0, s.t. if |A| < A i; £/ien 

dist(spec(-R A (a)), spec(e irL5 )) < C|A|r. (70) 

Moreover, the group of eigenvalues associated to any two distinct eigenvalues 
e lre 7 e lre of Rq{ol) = e lrLs belong to disjoint balls centered at e 1Te and e lre . 

A proof is obtained from a straightforward estimate of the resolvent 
(R\(a) — z)~ l using the Neumann series and the fact that || (Ro(a) — ^ ) 1 1 1 = 
[dist(z, speci?o(o ; ))]~ 1 since Ro(a) is normal. Due to (69), the motion of 
eigenvalues of R\(a) under variation of A is estimated for |A| < Ai, uniformly 
in a G C, |Imo;| < a , see [26] Section II §3. 
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